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Fuel prices are rising, and fossil fuel resources are depleting. This fuel is 
consumed by conventional fuel vehicles, which contributes to increased 
greenhouse gas emissions and air pollution. Various vehicle technologies 
have been introduced in modern times. Among these innovations, fuel cell 
vehicle technology stands out due to its ease of use, robustness, and simplicity. 
The primary energy source for this system is hydrogen. Fuel cells use 
hydrogen to generate electricity through a chemical process that does not 
involve combustion. The purpose of this paper is to describe the design and 
development of a simulation model for a fuel cell (proton exchange membrane 
type) vehicle using mathematical equations and the MATLAB/Simulink 
R2020a software, taking into account all voltage losses and temperature 
variations of the proton exchange membrane (PEM) type fuel cell. The 
designed model is validated in two ways: first by generating V-I 
characteristics for the designed PEM-type fuel cell, and then by performing a 
performance test in Simulink using the driving cycle on one of the fuel cell 
vehicles. The obtained simulated results are nearly identical to the expected 


results, and the designed PEM-type fuel cell performed well in all real-time 
driving cycle test conditions. 
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1. INTRODUCTION 

The Indian government is actively implementing stricter emission regulations to address the problem 
of vehicle pollution. Fuel price increases and the depletion of fossil fuel resources highlight the need for new 
technological advances in the automotive industry. Adoption of fuel cell hybrid vehicles (FCHVs) is one 
suggested solution to this problem. 

The burning of fossil fuels is the primary contributor to both localized air pollution and global 
warming. The transportation and automobile sectors are the largest consumers of fossil fuels; as a result, the 
primary objective is to reduce and eliminate pollution from the transportation and automobile sectors [1]. 
Vehicle manufacturers are competing against one another to produce environmentally friendly automobiles 
because this is a current necessity. Now that electric vehicles are available on the market, however, consumers 
should be aware that these automobiles have several drawbacks, including a short driving range, a lengthy 
amount of time required to recharge the battery, and inefficient operation. These drawbacks of electric vehicles 
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can be circumvented by using a fuel cell hybrid vehicle (FCHV), which operates on batteries and fuel cells that 
use hydrogen as a fuel [2], [3]. 

Researchers have conducted experiments on a wide variety of fuel cell designs. One of the potential 
candidates is a fuel cell of the proton exchange membrane (PEM) variety. The electrolyte that is used in a fuel 
cell is typically used to differentiate between the various types of fuel cells; however, there are always other 
significant factors to consider. The investigation of the mathematical modeling of PEM-type fuel cells using 
MATLAB/Simulink is beneficial because of its broad applicability, potential for efficiency improvement, and 
environmental advantages. When selecting a fuel cell technology, however, researchers and engineers should 
be aware of the cost and durability challenges associated with PEM-type fuel cells and should take into 
consideration the particular application and context [4], [5]. 

PEM-type fuel cells are most suitable for vehicle applications because they allow for rapid start-up 
thanks to their low operating temperature. Internal-combustion engines find these to be an attractive choice 
because of their high-power density and robust mechanical structure [1], [4], [6]. The objective of this study is 
to create and construct a simulation model for proton exchange membrane (PEM) fuel cells specifically 
designed for use in vehicles, particularly four-wheeled vehicles. The model is formulated using mathematical 
equations through the utilization of MATLAB/Simulink software. This software aids in testing and assessing 
the feasibility of the output power in accordance with the specified requirements. This will also aid in 
comprehending the distinct operational zones of PEM-type fuel cells and identifying the optimal zone for PEM- 
type fuel cells to function in relation to voltage and current. Creating a MATLAB/Simulink model for a proton 
exchange membrane (PEM) fuel cell is a flexible and highly valuable tool. It has versatile applications, such 
as enhancing performance, analyzing systems, developing controls, conducting research, and facilitating 
education in the domain of fuel cell technology. 


2. METHOD 
2.1. Proton exchange membrane (PEM) fuel cell 

A PEM-type fuel cell is made up of two electrodes: a negative-polarity anode and a positive-polarity 
cathode, as well as a membrane of solid electrolyte between them. It is the same as an electrochemical cell. 
Oxidization of hydrogen is on the anode side, while reduction of oxygen is on the cathode side. Anode-to- 
cathode proton transmission happens through the membrane of a solid electrolyte. Electrons flow to the cathode 
over an external circuit. On the cathode side, water and heat are generated as byproducts of a combined reaction 
between O2, electrons, and protons. The electro-chemical reaction is accelerated due to the presence of a 
catalytic member on the anode as well as the cathode [7], [8]. A typical PEM-type fuel cell has the following 
reactions [7]: 

— Anode: H2(g)—2H*(aq) +2e7 
— Cathode: %02(g) +2H*(aq) +2e —H20(1) 
— Overall: H2(g) +/202(g)—H20(1) + electrical energy + thermal energy (heat as a by-product) 

The fuel source for PEM fuel cells used in automotive applications is hydrogen gas. Figure 1 [9] 
shows how a fuel cell of the PEM type operates. The diagram shows that air is supplied to the cathode side and 
hydrogen gas passes through the anode side. The electrolyte in the fuel cell divides hydrogen into protons and 
electrons. Electrons follow an external circuit to create an electric current, while protons travel through the 
electrolyte to reach the cathode side. Heat and water are produced on the cathode side when oxygen, protons, 
and electrons mix. PEM-type fuel cells are suitable for automotive applications as they operate at low 
temperatures and give quick responses by varying output current and voltage to meet power requirements in 
vehicles as per load conditions. 


2.2. Mathematical modelling of PEM type fuel cell 
To model the fuel cell, the first step required is to calculate Nernst voltage of fuel cell and it is given 
by (1) [10]: 


AS RT 
Ey = Eq + (T — Ty) = — = In(PH, P02) (1) 


where, E,,= fuel cell’s Nernst voltage; Ey= standard-state reversible voltage; T = temperature of fuel cell in K; 
Tp = room temperature in K; AS = change in entropy for H2-O2 fuel cell; z = no. of electrons that were 
transported during the process; F = Faraday’s constant = 96485 Cmol"!; R= gas constant =8.314 J (molK) |; 
PH = partial pressure of hydrogen inside the stack; PO2 = partial pressure of oxygen inside the stack. 
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Figure 1. Working of a PEM type fuel cell [9] 


2.3. Identification of the Nernst voltage equation’s key elements 

A number of crucial factors need to be considered in order to calculate the Nernst voltage. Since 
hydrogen-oxygen fuel cells are our main area of interest, the following parameters are calculated specifically 
for this setting. A hydrogen—oxygen type fuel cell generates the standard-state reversible voltage (Eg) and is 
determined using (2) [10]. 


ee 
Ey=-F (2) 


Where Ag = standard-state free energy change for the hydrogen-oxygen reaction; Ag = —237000 J mol; z =2, 
for hydrogen-oxygen fuel cell. Substituting all the known values, in (2). 


Ey = —237000 
2 X96485 
Ey = +1.229V 


For an H»-O2 fuel cell operating in standard state conditions, the change in entropy (AS) is found 
to be —44.34 J (mol K)"!, where (Tp) is the room temperature, which is 298 K. The partial pressure of oxygen 
(PO2) and hydrogen (PH2) depend on the vehicle’s power requirement and reactant consumption. In the end, 
(3) represents the Nernst voltage equation for a hydrogen-oxygen fuel cell. 


E, = 1.229 — 2.298 x 10-* (T — 298) — = n(PH, /PO;) (3) 


2.4. Fuel cell voltage 
Voltage of one fuel cell and stack is determined by (4) and (5) respectively [11], [12]: 


Veet = Ey ~ Vact,cell a Vonmic.cell ~ Veonccell (4) 
Vetack = N X Veen (5) 


where, Ve; = Voltage o/p of one fuel cell; Vact.cey = activation loss of one fuel cell; Vopmiccen = Ohmic loss of 
one fuel cell; Veonecen = concentration loss of single fuel cell; Vetack = fuel cell stack voltage; N = Number of 
cells in stack. The number of cells, required to obtain the anticipated output power, is computed using (6). 


Voltage output required 


i Voltage output of single fuel cell (6) 
2.5. Voltage losses in PEM type fuel cell 

Due to irreversible losses, the present fuel cell’s voltage output is observed as lower than that 
estimated by thermodynamics principles. As shown in Figure 2, the following are three significant forms of 
fuel cell losses: in early part of curve, the electro-chemical process causes activation losses. 
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Ideal (thermodynamic) fuel cell voltage 
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Figure 2. PEM type fuel cell voltage variation with net current density [11] 


The ionic electronic state causes ohmic losses, in the central part of curve. Mass transfer causes concentration 

losses in the end of V-I curve [11]. 

a) Activation loss: the phenomenon known as activation loss is the slow rate at which reactions take place on 
the electrode surfaces. As a result, some of the voltage produced is highly nonlinearly dissipated in order 
to support the chemical reactions that carry out the electron transfer to and from the electrode. The 
activation loss of fuel cell is given by (7) [4], [13]. The ‘Tafel equation’ is the name given to (7). The (8) 
helps to determine the Tafel slope. 


Vact = AX In(>) (7) 
mes 
eT (8) 


Where, A = Tafel slope (V); I = current density (Ampere cm”); i, = exchange current density 
(ampere cm~); a = charge transfer coefficient in the range 0 to 1.0. 

Ohmic loss: the electrodes' electrical resistance is the cause of ohmic losses. The movement of ions in the 
electrolyte is hindered by a variety of connections and resistances. There is a linear relationship between 
the current density and the magnitude of the voltage drop. The ohmic loss of fuel cell is 
given by (9) [4], [12]. 


b 


wm 


Vonm = 1 X Rohm =! X (Relectronic + Rionic) (9) 


Where, / = fuel cell current (ampere); Ronm = ohmic resistance (ohm ‘Q’); Rejectronic = resistance to the 
electron transport in electron conductive part (ohm ‘Q’); Rionic = resistance to the flow of ions (ohm ‘Q”). 
c) Concentration loss: concentration gradients are created during the reaction process as a result of the mass 
diffusion of reactants from the flow channels to the reaction sites, which causes concentration loss. The slow 
transport of reactants to the reaction sites causes concentration voltages to drop at high current densities, 
resulting in a non-linear voltage drop. The concentration loss of fuel cell is given by (10) [4], [14]. 


RT I 
Voone = oF In(1 — i (10) 
Where J, = limiting current of fuel cell (ampere). 


2.6. Thermodynamic energy balance 
The net heat generated inside the PEM type fuel cell due to chemical reaction causes its temperature 
to increase or decrease. This net heat generation given by (11) [15]. 


Anet = AUchem— elec — Ysens+latent — Moss (11) 


Where, dnet = net heat generated inside the PEM type fuel cell; qchem = chemical energy generated inside the 
PEM type fuel cell; qdelec = electrical energy generated inside the fuel cell; dsens+jatent = the sensible and 
latent heat generated inside the fuel cell; q}9;, = heat loss due to convection inside the fuel cell. 
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2.7. Designed PEM type fuel cell 

A PEM-type fuel cell is designed, modeled, and simulated using MATLAB/Simulink software. The 
designed fuel cell is modeled with the help of mathematical equations. The advantage and uniqueness of this 
designed PEM-type fuel cell model is that it consists of all PEM-type fuel cell voltage losses along with 
thermodynamic energy balance. As shown in Figure 3, the MATLAB model has mainly three subsystems. 
Nernst voltage and rate of utilization are calculated in one subsystem. In another second subsystem, all voltage 
losses are calculated for a single fuel cell, and then losses across the fuel cell stack are calculated by considering 
the number of cells. In the third thermodynamic energy balance subsystem, the temperature of the fuel cell is 
calculated. Voltage losses block as depicted in more detail in Figure 4 gives the activation voltage drop, ohmic 
voltage drop, and concentration voltage drop of a single cell. Figure 5 depicts the detailed simulation block 
pertaining to determination of rate of utilization of reactant and oxidant i.e., H2, O2 and O2H2. Other inputs 
required to calculate Nernst voltage are anode and cathode side supply pressure, and the composition of 
hydrogen and oxygen. Figure 6 depicts that Tafel slope, exchange current density and change in enthalpy are 
also calculated along with Nernst voltage. Figure 7 depicts the thermodynamic energy balance subsystem. This 
subsystem helps to calculate net heat generated inside the fuel cell, and further helps to determine temperature 
of fuel cell. 
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Figure 4. Subsystem voltage losses block 
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Figure 5. Subsystem rate of utilization block 
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Figure 7. Subsystem thermodynamic energy balance 
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3. RESULTS AND DISCUSSION 
3.1. Validation of designed PEM type fuel cell 

The objective of the research is to carry out mathematical modelling of PEM type fuel cell for vehicle 
applications (4-wheeler) which further helps to test and get the idea of output power feasibility as per 
requirement. To design PEM type fuel cell, required parameters as mentioned in earlier sections are decided 
as: 1) The fuel cell must produce 100 kW of power; ii) The fuel cell's operating temperature range is 0 to 
100 °C; iii) The cathode and anode sides are under 3 atmospheres of pressure, or 303,975 Pascals (Pa); and 
iv) Number of cells calculated by using in (6) are 400. Based on the previously mentioned parameters and the 
equations presented in the earlier sections, a MATLAB/Simulink mathematical model of a PEM-type fuel cell 
is built following the methodology described. 

Similar to the V-I characteristics shown in Figure 2, Figure 8 shows the voltage-current (V-I) 
characteristics of the PEM-type fuel cell. Using a MATLAB/Simulink software-available representative fuel 
cell vehicle example, the PEM-type fuel cell is tested and validated [16], [17]. The PEM-type fuel cell has a 
maximum power demand of 100 kW, and the maximum power output, according to the data, is also 100 kW. 
With a decreasing voltage as the current increases, the PEM-type fuel cell's performance closely resembles the 
fuel cell's theoretical behavior. The voltage of the fuel cell operates in three different regions: concentration, 
ohmic, and activation. The induced voltage in the PEM-type fuel cell that has been designed ranges from 481.2 
to roughly 390 volts during the activation phase. The slow chemical reaction within the activation region causes 
a voltage reduction, resulting in a non-linear voltage drop. 

Immediately after the activation region ends, the fuel cell begins to function in the ohmic region. This 
voltage ranges from about 390 volts at the beginning to about 315 volts at the end. Electrical resistance in the 
electrode causes ohmic loss within this ohmic region, which results in a linear voltage drop. It's important to 
note that this area, which corresponds to the point of maximum performance at 210 amperes and 340 volts, is 
thought to be ideal for fuel cell operation. The fuel cell then moves from the ohmic region to the mass transport 
region, which begins at about 310 volts and lasts until the cycle ends, when it approaches 290 volts. Diffusion 
of reactants from the flow channels to the reaction sites causes mass concentration loss in the mass transport 
region. Because reactants are transported to the reaction sites more slowly in higher-current operating zones, 
this loss is more noticeable. 
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Figure 8. Simulated VI characteristics of the designed PEM type fuel cell 


3.2. Performance test of designed PEM type fuel cell using NEDC 

The performance test of the designed PEM-type fuel cell is done with the help of ‘New European 
Driving Cycle’ (NEDC) [18]. The NEDC is given as input to one of the examples of fuel cell vehicles available 
in MATLAB/Simulink software. NEDC provides different test points at different speeds, and the power 
requirement and performance test of the vehicle for different real-time test conditions are observed. 

The NEDC cycle runs for 1160 seconds. The output voltage, output current, efficiency, temperature, 
and output power performance of the designed PEM-type fuel cell are observed. During performance testing, 
the output voltage and current behavior [19] of the designed PEM-type fuel cell follow those of the theoretical 
fuel cell. As the speed of the vehicle increases, the current demand increases, and as the current increases, the 
fuel cell’s voltage decreases. After 800 seconds, the vehicle's speed increased, and further speed increased to 
120 kmph, as illustrated in Figure 9. The output current of the designed PEM-type fuel cell, as illustrated in 
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Figure 10(a), is increasing, and in that period, the output voltage, as shown in Figure 10(b), is decreasing to 
satisfy the power requirement of the vehicle through the motor of the vehicle. 

The fuel cell's irreversible efficiency is in the range of 40 to 60% [20]. The efficiency of the designed 
PEM-type fuel cell can be observed in Figure 11(a). By ignoring the lowest and highest peaks, the mean 
efficiency of a PEM-type fuel cell, as observed in Figure 11(a), is also in the range of 45 to 60%. The fuel 
cell’s temperature operating range is from 25 to 100, i.e., 298 K to 373 K. As observed from Figure 11(b), the 
range of temperatures within which the designed PEM-type fuel cell is operating is also 298 to 373 K, which 
is near the stated operating temperature range of the fuel cell. 

The vehicle's electric motor efficiently supplies the required power depending on a number of 
variables, including speed and load. The vehicle's battery and fuel cell power work together to meet these 
power needs. The PEM-type fuel cell intended for the NEDC cycle is tested, and its performance is assessed 
under a range of real-time test conditions in Figure 12. The designed PEM-type fuel cell and the vehicle's 
battery work together to guarantee that the electric motor's power needs are sufficiently satisfied. As observed 
from Figure 12, the motor power requirement of the motor is up to 25 kW, supplied by the designed PEM type 
fuel cell, and the remaining power of the designed PEM type fuel cell is used for regenerative braking, i.e., for 
the charging of the battery when it is not supplying power. The motor power requirement at high speed is from 
1000 seconds to 1160 seconds, i.e., from 30 kW to about 70 kW. Thus, when the designed PEM-type fuel cell 
alone cannot fulfill the power requirement of the motor, the remaining power requirement of the motor is 
supplied by the battery of the vehicle, which is charged during regeneration mode. In this way, the validation 
of the designed PEM-type fuel cell is done using NEDC. 
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Figure 9. New European driving cycle (NEDC) 
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Figure 10. Performance of designed PEM type fuel cell (a) output current and (b) output voltage 
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Figure 11. Designed PEM type fuel cell, (a) efficiency and (b) temperature 
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Figure 12. Output power at different load and speed conditions by designed PEM fuel cell and 
battery for motor 


Fuel cell modeling can be done by considering the double-layer charge effect [21]. In the double-layer 
charging effect, at the junction point where the porous cathode meets the membrane, two polar, oppositely 
charged layers are formed. These are electrochemical double layers [4], [22]. The fuel cell model incorporates 
the double-layer charged effect by calculating output voltage using concentration voltage drop [23]. 

When the frequencies are higher, the impedance in a fuel cell acts like a resistance with an inductance 
that is connected in series. Because the double-layer capacitor filters current harmonics at these frequencies, 
they have negligible influence in regard to the reactions that are chemical in nature [24], [25]. The lifetime of 
a PEM-type fuel cell is influenced by low-frequency current ripple, which must be controlled within a certain 
minimum range. To control the low-frequency current ripple, several control methods, such as current mode 
control and current mode control with a notch filter in the voltage loop, are necessary [26], [27]. As the scope 
of this research is limited to the modeling of the PEM-type fuel cell, further work can be done to reduce low- 
frequency current ripple by implementing various control methods. 


4. CONCLUSION 

Mathematical modeling of PEM-type fuel cells is done for vehicle applications (4-wheelers). Fuel cell 
models are available in MATLAB and Simulink, and by using these examples, this PEM-type fuel cell model 
is validated. The results drawn from the designed PEM-type fuel cell model are as follows: i) This model of a 
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PEM-type fuel cell has the same V-I (voltage-current) characteristics as the PEM-type fuel cell that is thought 
to work in theory; ii) A PEM-type fuel cell meets the maximum power requirement. Furthermore, at various 
load conditions, the proposed PEM-type fuel cell's power output is in accordance with the vehicle's power 
requirements; iii) The best operating region for PEM-type fuel cells is the ohmic range, where PEM-type fuel 
cells operate within their optimum range with better efficiency; iv) NEDC was used to test the behavior of the 
designed PEM-type fuel cell. The stack properties observed are comparable to those achieved in prior 
investigations. The behavioral statistics related to energy consumption are the same as those from earlier 
studies; v) The designed PEM-type fuel cell’s efficiency is between 45% and about 60%; vi) The temperature 
of the designed PEM-type fuel cell is within the range of 25 to 100, which is as per the acceptable requirement 
for functioning. 
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